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A nonlinear screening theory is formulated to study the problem of gap formation and its relation 
to glassy freezing in classical Coulomb glasses. We find that a pseudo-gap ("plasma dip") in a 
single-particle density of states begins to open already at temperatures comparable to the Coulomb 
energy. This phenomenon is shown to reflect the emergence of short range correlations in a liquid 
(plasma) phase, a process which occurs even in the absence of disorder. Glassy ordering emerges 
when disorder is present, but this occurs only at temperatures roughly an order of magnitude lower. 
Our result demonstrate that the formation of the "plasma dip" at high temperatures is a process 
distinct from the formation of the Efros-Shklovskii (ES) pseudo-gap, which in our model emerges 
only within the glassy phase. 



The interplay of interactions and disorder remains one 
of the most important open problems in condensed mat- 
ter physics. These effects are most dramatic in disor- 
dered insulators, where the pioneering work of Efros and 
Shklovskii (ES) [l| emphasized the fundamental signifi- 
cance of the long-ranged nature of Coulomb interactions. 
This work presented convincing evidence that at T = a 
soft " Coulomb gap" emerges in the single-particle density 
of states (DOS) which, in arbitrary dimension d, reads 
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^From a general point of view this result is quite sur- 
prising. It indicates a power-law distribution of exci- 
tation energies, i.e. the absence of a characteristic en- 
ergy scale for excitations above the ground state. Such 
behavior is common in models with broken continuous 
symmetry, where it reflects the corresponding Goldstone 
modes, but is generally not expected in discrete symme- 
try models, such as the one used by ES. Here, it may 
reflect unusually strong frustration behavior inherent to 
Coulomb interactions in presence of disorder. 

Indeed, the ES model seems to display several glassy 
features characterized by a large number of meta-stable 
states and slow relaxation, as clearly seen in many simu- 
lations y, y, y, y| , and even in some experiments jZi] . In- 
terestingly, a precursor of the gap begins to appear [J, |g 
already at relatively high temperatures, while glassy or- 
dering emerges only at much lower T. Similar behavior 
has been identifled even in absence of randomness |g . Is 
the physics of the Coulomb gap thus related or is it unre- 
lated to the glassy features of the system? The close con- 
nection between the two phenomena was recently demon- 
strated loj for a mean-fleld model of interacting disor- 
dered electrons in the limit of large coordination, but the 
issue remains unresolved for Coulomb systems in finite 
dimensions where the ES theory applies. 



To address these important issues in a controlled and 
precise fashion, we make the following observation which 
is the main physical point of this letter. We stress that 
the principal ES result - the emergence of a power-law 
spectrum - is not specific to low dimensions! Its phys- 
ical origin and its relation to high temperature anoma- 
lies can, therefore, be investigated by the theoretical ap- 
proaches controlled in the limit of high spatial dimen- 
sions. Such a theory for the Coulomb glass (CG) model 
is presented in this letter. Our main conclusions are as 
follows: (1) A non-universal pseudo-gap in the DOS (we 
call it the "plasma dip") begins to emerge at temper- 
atures of the order of the Coulomb energy. It reflects 
strong short range correlations in the Coulomb plasma, 
a feature that is most pronounced in absence of disorder 
[3, but is unrelated to the Coulomb gap of ES. We ob- 
tain simple analytical results that in quantitative detail 
describe the temperature evolution of the plasma dip, 
in excellent agreement with all existing simulations. (2) 
The high temperature plasma (fluid) phase becomes un- 
stable to ergodicity breaking at temperatures typically 
ten times lower, as the system enters a glassy state. We 
argue that a true ES pseudo-gap emerges only within the 
glassy (nonergodic) phase, and that its scale-invariant 
form reflects the marginal stability of such a glassy state. 

Nonlinear screening theory. The simplest many-body 
approach to Coulomb systems is the well-known Debye- 
Huckel theory , which provides a linear screening descrip- 
tion equivalent to a Gaussian approximation for the plas- 
mon mode. This formulation, however, fails badly at low 
temperatures, where nonlinear effects lead to strong cor- 
relations in the plasma phase. Even worse, such a Gaus- 
sian theory is unable to describe glassy freezing even in 
the well-understood mean-field limit corresponding to in- 
finite range interactions. To overcome these difficulties, 
we use the simplest theory of nonlinear screening given 
by the classical limit of the so-called extended dynamical 



mean-field theory [lOj , which also describes the leading 
order nontrivial correlations in the limit of large coor- 
dination. In this approach, the environment of a given 
site is approximated by free collective modes (plasmons 
in our case), the dispersion of which is self-consistently 
determined. In recent work, a version of this method has 
been successfully applied to the problem of self-generated 
glassiness [ll| in systems with frustrated phase separa- 
tion [lj| without disorder. Here, we apply it to the clas- 
sical Coulomb glass model ^ij given by the Hamiltonian 



H = Y, 4>^n^ + IY1 ^^1 ("* - ^)("J - ^)' 
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where rii = 0,1 is the electron occupation number, and 
(j)i is a Gaussian distributed random potential of variance 
W^. We express the Coulomb interaction Vij — ea/rij in 
units of the nearest- neighbor repulsion eoi a-nd the inter- 
site distance r^ in units of the lattice spacing. We adopt 
the following notation throughout the paper. For vectors 
and matrices in the replica space we use bold font and 
the hat symbol correspondingly. In the replica symmetric 
ansatz (RS) a matrix O = {Oc, 0} is parameterized by its 
connected part Oc and its off-diagonal part O. Thus for 
the density density correlator we use g = {Xi 9}- Thermal 
and disorder averages are denoted as {0)t and [0]dis 
respectively, and (O) = [{0)T]dis- 

To derive the desired self-consistency equations we av- 
erage over disorder using the standard replica method [2| , 
and use a cavity construction |2|, |l(j| , integrating out the 
degrees of freedom on all sites except the considered one. 
The resulting contribution to the local effective action is 
computed in the Gaussian approximation, giving a term 
of the form — -icJnAtJn, where 5n — n— (n). By requiring 
that the local density-density correlator q°'P = {6n°'5n^)c 
is correctly reproduced by the effective action, one ob- 
tains a self-consistency condition. We get 



g"/3 = {5n°'5n'^) 
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where A"'' = A"'' + P^W^, and Vk is the Fourier trans- 
form of the interaction potential. In case of the RS solu- 
tion A 
the renormalized disorder. 



{A,, (3^W^ff}, where Weff 



y^W^ + /3-2A is 



Fluid solution. To examine the evolution of the system 
in the high temperature phase we first examine the RS 
solution. In the n — > replica limit we get 
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where D[x] = (2Tr)-^/'^ exp{-x'^/2} dx. The self- 
consistency condition becomes 
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The Eqs. H4I5(I can easily be solved numerically for 
the parameters XiQT^c,A,Weff, to calculate the den- 
sity of states (DOS) function g{e), to examine the sta- 
bility of the RS solution, and to compute the entropy 
S. The above equations are written for the half filled 
case in absence of charge ordering. Generalization to uni- 
formly ordered phases is a standard procedure [l^j , where 
the ordering transition is signaled by a divergence in 
Xk = {x~^ + Ac — /?Vfe) at the ordering vector k = Q. 
For simplicity, most of our results are written in the ho- 
mogeneous phase at half filling and that will be assumed 
unless stated otherwise. 

Density of states. For the CG model, the single particle 
DOS (tunneling DOS) function g{e) is simply given by 
the distribution of the local fields (energies) e^ = dH/dni. 
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Integrating out all sites except one, we derive an expres- 
sion for g{e) in terms of the local effective field A. The 
final result for the RS solution reads: 



g{e) = / D[x 
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In the cavity method langauge 14^ this result can be in- 
terpreted as the Gaussian distribution of the cavity fields 
(different from the local fields) of variance We// , modi- 
fied by the the (Onsager) self reaction term representing 
the plasma correlations. 

We get insight into the behavior of the DOS by con- 
sidering some analytically solvable limits. At T — > 0, the 
DOS remains finite at the Fermi level, though it can be 
exponentially small for weak disorder. This shows that 
the fluid (RS) solution does not capture the physics of 
the true ES gap, which emerges only within the glassy 
phase. On the other hand, as long as the RS solution is 
stable, a large "plasma dip" may develop, but it will have 
no direct relation to the glassy physics or the ES gap. It 
reflects strong short-range correlations in the Coulomb 
plasma (fluid) phase, which are suppressed at large dis- 
order. Here, the RS DOS reduces to the bare disorder 



distribution. In the opposite limit of vanishing disorder 
Weff -^ the DOS expression simphfies and for an ar- 
bitrary fiUing reads: 
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To support the vahdity of our formulation, we com- 
pare our analytical results with available numerical sim- 
ulations. In Fig. ^ (top panel), we examine the situation 
studied in Ref . |a] , where calculations were done for a 3D 
CG on a cubic lattice, for a set of temperatures in the 
fluid phase. The lowest temperature is very near the glass 
transition temperature (see our phase diagram. Fig. 2). 
As we can see, our theory captures in surprisingly quan- 
titative detail the formation of the "plasma dip" in the 
fluid phase. In the past, this phenomenon has often been 
confused with the formation of the true ES gap which, 
as we argue below, only emerges within the glassy phase. 
Finally, we test limits of our theory by computing the 
DOS for the 2D CG in absence of disorder. Even in this 
extreme case, we reproduce semi-quantitatively exact nu- 
merical results of Ref. |3| • 

Glassy ordering. To examine the stability of the fluid 
phase to glassy ordering, we examine the Baym Kadanoff 
(BK) functional Tbk [q] ■ This is a functional of the cor- 
relator q, which yields the exact equations of motion at 
the saddle point, where it coincides with the exact free 
energy. To obtain our self-consistency conditions, a local 
approximation 13] is made on the two particle irreducible 
part oiTsK- In this formulation, the stability of our fluid 
RS solution can be obtained by a standard replica sym- 
metry breaking (RSB) analysis T^j of the BK functional 
at the saddle point. The corresponding RSB instability 
criterion takes the form 
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Here, Xij is the density-density correlation function com- 
puted for a fixed realization of disorder, i.e. Xij — 
{Snidnj)T — {Sni)T{Snj)T- The left hand side of Eq. ^ 
is nothing but ^/J2j[Xij]dis, the inverse of the glass sus- 
ceptibility, a quantity which diverges at the transition. 
In terms of the RS solution the RSB condition reads 
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As an illustration, we present results for the CG on a 3D 
cubic lattice, and in Fig. |21 we plot the corresponding 
phase diagrams obtained by numerically solving our self- 
consistency conditions. At small disorder and tempera- 
ture T w .95 (which is in satisfactory agreement with the 
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FIG. 1: Our analytical predictions for the single-particle den- 
sity of states (full lines) are found to be in excellent quanti- 
tative agreement with simulation results (dashed lines), with 
no adjustable parameters. Shown are results for the three 
dimensional case studied in from Ref. 1!], corresponding to 
W = 1/(273), and temperatures T = 0.4,0.2,0.1,0.05 (top 
panel), and the two dimensional model of Ref. [g, corre- 
sponding to VF = 0, T = 0.1, K = 0.2. All lines correspond 
to the plasma phase, while T — 0.05 is very close to the glass 
transition temperature. 



exact value|l7j Tc = .129) the system enters the charge 
ordered phase, phase. Stronger disorder suppresses the 
charge ordering, and the system can exist either in a liq- 
uid phase (at higher temperature) or in the glass phase 
(at lower T). The liquid is separated from the glass by 
the RSB line, also known as the Almeida-Thouless [l5j 
line. We emphasize that the ordering temperature we 
predict is roughly an order of magnitude smaller then 
the Coulomb energy, in remarkable quantitative agree- 
ment with all available simulation results |^ Ij, M, |^ . 

This interesting fact can be traced down to the screen- 
ing of the Coulomb interaction. Indeed, the overall en- 
ergy scale characterizing the screened Coulomb potential 
Vscrir) — Eq exp{—r/£scr}/r is roughly an order of mag- 
nitude smaller then the bare Coulomb energy. The cor- 
responding screening length €scr = [(x~^+Ac)/(/3eo)]^''^ 
(shown for W = (2V5)~^ in the inset of Fig. 2) de- 
creases (albeit weakly) with temperature, and remains 
short throughout the fiuid (RS) phase. This observation 
also makes it clear why the true ES gap cannot emerge 
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FIG. 2: 3D Coulomb glass phase diagram. The full horizontal 
line indicates the RSB instability and the dotted line shows 
where the RS entropy turns negative. The screening length 
in the inset is plotted for the same disorder and the range of 
temperatures (fluid phase), as in Fig0 



in absence of glassy ordering. The screening mechanism 
remains operative throughout the (ergodic) fluid phase, 
and thus the long-range character of the Coulomb in- 
teraction, which is crucial for the ES argument remains 
inoperative. In contrast, within the glass phase, follow- 
ing the arguments from Ref. '9f|, we expect the relevant 
zero-field cooled compressibility to decrease and vanish 
at T = 0. This mechanism opens a route for the screen- 
ing to be suppressed at low temperatures, and the true 
ES gap to emerge. 

To provide further evidence of the instability of the 
fluid phase to glassy ordering, we also calculate the en- 
tropy in the fluid (RS) solution, which takes the form 



S= / D[x]ln 



2 cosh ( ^x(3W,ff 



+ ^J2^nxk-llnx-2xl3W,ff. (11) 



It is well known that for standard mean-field glass mod- 
els, the RS replica theory predicts negative entropy at 
T = 0, while the lower bound for the glass transition 
temperature can be set where the RS entropy changes 
sign. It is easy to show from Eq. (|ll|l that our RS en- 
tropy proves strictly negative at T = as well. In Fig. 
2 we also plot (dashed line) a lower bound for the RSB 
instability line where the RS entropy turns negative, pro- 
viding further evidence that the fluid phase cannot sur- 
vive down to T = 0. 

The glass phase and the Efros-Shklovskii gap. In this 
letter we do not explicitly examine the RSB solution of 
our model. Nevertheless, we follow arguments similar to 
those of Ref. [9| , and use the large disorder asymptotics 
of the glass transition line to determine the powerlaw 



form of the T = Coulomb gap in the glassy phase. At 
W ^ eo, we find 



Tn 



e^V- 



(12) 



where the exponent 7 is identical to that predicted by 
ES for the T = DOS in the CG 
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For an interaction of a general power law form V(r) ~ 
Sa/r"", the ES argument predicts 7 = (d — a)/a in arbi- 
trary dimension d. The asymptotic regime, however, sets 
in at larger values of disorder, not shown in the Fig.©. 
Let us explain the importance of Eq. H12|) . We have 
seen that for W y^ 0, the DOS remains finite in the RS 
fluid phase, so a true ES pseudo-gap can emerge only 
due to glassy ordering. For another electron glass model, 
work of Ref. [2j has established that the emergence of 
a true pseudo-gap at T = directly follows from the 
marginal stability of the glassy state, and that its form 
also determines the large disorder asymptotics of Tg(W). 
Given the close similarity of our mean-field equations for 
the CG model to those examined in Ref. .^], we expect 
the same mechanism to apply here as well. Using the 
expected form of the ES gap, we can estimate the glass 
transition temperature as the energy scale Egap charac- 
terizing the "width" of the gap that opens in the low 
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temperature phase. We find E, 



gap 



w- 



coincid- 



ing with Eq. (|12|l . This result presents strong evidence 
in favor of the close relation between glassy ordering and 
the emergence of the ES gap. 

In conclusion, we have formulated a simple many-body 
theory that is able to clarify the relation between the fi- 
nite temperature formation of the Coulomb gap and the 
emergence of glassy ordering in disordered Coulomb sys- 
tems in finite dimensions. This nonlinear screening ap- 
proach is fiexible enough to allow for future extensions 
to quantum models ^ J^] , and to study the role of An- 
derson and Mott localization [l3| in Coulomb systems. 
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